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Abstract 

It is shown that prize changes of the US dollar - German Mark ex- 
change rates upon different delay times can be regarded as a stochastic 
Marcovian process. Furthermore we show that from the empirical data 
the Kramers-Moyal coefficients can be estimated. Finally, we present an 
explicite Fokker-Planck equation which models very precisely the empirical 
probabilitiy distributions. 

PACS: 02.50-r;05.10G 

Since high-frequency intra-day data are available and easy to access, research 
on the dynamics of financial markets is enjoying a broad interest. ]1|, [2|, ||, [4], 
H> B 01- Well-founded quantitative investigations now seem to be feasible. The 
identification of the underlying process leading to heavy tailed probability density 
function (pdf) of price changes and the volatility clustering (see fig.l) are of special 
interest. The shape of the pdf expresses an unexpected high probability of large 
price changes on short time scales which is of utmost importance for risk analysis. 
In a recent work [§], an analogy between the short-time dynamics of foreign 
exchange (FX) market and hydrodynamic turbulence has been proposed. This 
analogy postulates the existence of hierarchical features like a cascade process 
from large to small time scales in the dynamics of prices, similar to the energy 
cascade in turbulence c.f. 0. This postulate has been supported by some work 
M and questioned by others M . 
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One main claim of the hypothesis of cascade processes is that the statistics of 
the time series of the financial market can be determined. The aim of the present 
paper is to discuss a new kind of analysis capable to derive the underlying math- 
ematical model directly from the given data. This method yields an estimation 
of an effective stochastic equation in the form of a Fokker-Planck equation (also 
known as Kolmogorov equation). The solutions of this equation yields the prob- 
ability distributions with sufficient accuracy (see fig. 1). This means that our 
method is not based on the conventional phenomenological comparison between 
models and several stochastic aspects of financial data. Our approach demon- 
strates how multiplicative noise and deterministic forces interact leading to the 



heavy tailed statistics. [10[ 



Recently it has been demonstrated that from experimental data of turbulent 
flow the hierarchical features induced by the energy cascade can be extracted in 
form of a Fokker-Planck equation for the length scale dependence of velocity fluc- 



tuations [12[ . In the following it will be shown the dynamics of foreign exchange 
rates can also be described by a Fokker-Planck equation. Here we use a data set 
consisting of 1,472,241 quotes for US dollar-German mark exchange rates from 
the years 1992 and 1993 as used in reference f§. 

We investigate the statistical dependence of prize changes Axi := x(t + Ati) — 
x(t) upon the delay time At*. Here x(t) denotes the exchange rate at time t. We 
will show that the prize changes Axi, Ax 2 for two delay times Ati, Ai 2 are 
statistically dependent, provided the difference Ai 2 — Ati is not too large and 
that the shorter time interval is nested inside the longer one. 

In order to characterize the statistical dependency of price changes Axi, we 
have evaluated the joint probability density functions 

p(Ax 2 ,A* 2 ;Axi,A*i) (1) 

for various time delays At 2 < Ati directly from the given data set. One example 
of a contourplot of the logarithms of these functions is exhibited in figure 2. If the 
two price changes Axi, Ax 2 were statistically independent, the joint pdf would 
factorize into a product of two probability density functions: 

p(Ax 2) At 2 - Axi, Atx) = p(Ax 2) At 2 )p(Ax u Ati). (2) 

The tilted form of the probability density p{Ax 2) Ai 2 ; Axi, At\) clearly shows 
that such a factorization does not hold and that the two price changes are 
statistically dependent. This dependency is in accordance with observations 
of Miiller et al. for cross-correlation functions of the same data 0. To an- 
alyze these correlations in more detail, we address the question: What kind 
of statistical process underlies the price changes over a series of nested time 
delays Ati of decreasing length? In general, a complete characterization of 
the statistical properties of the data set requires the evaluation of joint pdfs 
p N (Axi, Ati] Axat, Ai^y) depending on TV variables (for arbitrarily large N). 
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In the case of a Markov process (a process without memory), an important simpli- 
fication arises: The N-point pdf p N are generated by a product of the conditional 
probabilities p(A%i, At^ + i|Aa^, At^), for i — 1, N — 1. As a necessary condi- 



tion, the Chapman-Kolmogorov equation [13 



p(Ax 2 ,At 2 \Ax u At 1 ) = J d(Ax i )p(Ax2,At 2 \Ax i ,At i )p(Ax i ,At i \Ax 1 ,At 1 ) 

(3) 

should hold for any value of Ati embedded in the interval 

At 2 < At 3 < At x . (4) 

We checked the validity of the Chapman-Kolmogorov equation for different 
Ati triplets by comparing the directly evaluated conditional probability distribu- 
tions p(Ax 2 , At 2 |Axi, Ati) with the ones calculated (p ca i) according to (|3|). In 
figure 3, the contour lines of the two corresponding pdfs are superimposed for 
the purpose of illustration; the red lines corresponding to p ca \. Only in the outer 
regions, there are visible deviations probably resulting from a finite resolution of 
the statistics. 

As it is well-known, the Chapman-Kolmogorov [|T3|] equation yields an evo- 
lution equation for the change of the distribution function p(Ax, At) across the 
scales At. 

For the following it is convenient (and without loss of generality) to consider 
a logarithmic time scale 

r = -InAt . (5) 

Then, the limiting case At —> corresponds to r —> oc. The Chapman- 
Kolmogorov equation formulated in differential form yields a master equation, 
which can take the form of a Fokker-Planck equation (for a detailed discussion 
we refer the reader to |l3fl): 



^-p(Ax,r) = 1-J^ D U(Ax,t) + J^DW(Ax,t)] P (Ax,t) (6) 
ar 2 oAx oAx z 

The drift and diffusion coefficients D^\Ax^ r), (Ax , r) can be estimated di- 
rectly from the data as moments of the conditional probability distributions 
(c.ffig. 3): 

dW(Ax,t) = him AT ^ M^ (7) 

M(k) = -^J dAx '( Ax> - Ax)V(Ax',r + At\Ax,t). 

As indicated by the functional Ax dependency of the moments , it turns out 
that the drift term is a linear function of Ax, whereas the diffusion term 
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D (2) 

is a function quadratic in Ax. In fact, from a careful analysis of the data 
set we obtain the following approximation: 

D m = -O.UAx 

D {2) = 0.003 exp (-r/2) + 0.019(Ax + 0.04) 2 (8) 

(Ax is measured in units of the standard deviation of Ax at At = 40960s). 
With these coefficients we can solve the Fokker-Planck equation (Q) for the pdf 
at times r > r with a given distribution at r . Figure 1 shows that the solutions 



of our model nicely fit the experimentally determined pdfs [|T4|. In contrast to 
the use of phenomenological supposed fitting functions c.f || ||, we obtained the 
changing forms of the pdfs by a differential equation. This method provides the 
evolution of pdfs from large time delays to smaller ones. This definitely is a new 
quality in describing the hierarchical structure of such data sets. At last, it is 
important to note that our finding of the Fokker-Planck equation for the cascade 
is in good agreement with the previously found phenomenological description in 



H . Based on the given exact solution of our Fokker-Planck equation [|TT| , we see 
that the chosen type of fitting function for the pdfs in || was the correct one. 
Furthermore, we see that the observed quadratic dependency of the diffusion term 
D 2 corresponds to the found logarithmic scaling of the intermittency parameter 
in ||, which was taken as an essential point to propose the analogy between 
turbulence and the financial market. 

We remind the reader that the Fokker-Planck equation is equivalent to a 
Langevin equation of the form (we use the Ito interpretation |T3H): 



^Ax(t) = DW(Ax(r),r) + Jd&)(Ax(t),t)F(t) . (9) 
dr v 

Here, F{r) is a fluctuating force with gaussian statistics 5-correlated in r: 

< F(t) >= , < F{r)F{r f ) >= 5(r - r f ) (10) 

In our approximation (§) the stochastic process underlying the prize changes is 
a linear stochastic process with multiplicative noise, at least for large values of 
Ax: 

^-Ax(t) = -0.44Ax(r) + v^019Ax (r)F(r) . (11) 
dr 

This stochastic equation yields realizations of prize changes Ax(r), whose ensem- 
ble averages can be described by the probability distributions p(Ax, r). Thus, the 
Langevin equation (§) produces the possibility to simulate the price cascades for 
time delays from about a day down to several minutes. Furthermore, with this 
last presentation of our results it becomes clear that we are able to separate the 
deterministic and the noisy influence on the hierarchical structure of the finance 
data in terms of the coefficients and D^ 2 \ respectively. 
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Summarizing, it is the concept of a cascade in time hierarchy that allowed us to 
derive the results of the present paper, which in turn quantitatively supports the 
initial concept of an analogy between turbulence and financial data. Furthermore, 
we have shown that the smooth evolution of the pdfs down along the cascade 
towards smaller time delays is caused by a Markov process with multiplicative 
noise. 

Helpful discussions and the careful reading of our manuscript by Wolfgang 
Breymann, Shoaleh Ghashghai and Peter Talkner are acknowledged. The FX 
data set has been provided by Olsen & Associates (Zurich). 
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figure captions 

Figure 1: Probability densities (pdf) p(Ax ) At) of the price changes Ax = 
x(t + At) - x(t) for the time delays At = 5120, 10240, 20480, 409605 (from bot- 
tom to top). Symbols: the results obtains from the analysis of middle prices of 
1,472,241 bit-ask quotes for the US dollar-German Mark exchange rates from 1 
October until 30 September 1993. Full lines: results form a numerical iteration of 
an effective Fokker-Planck equation with the initial condition of the probability 
distribution for At = 40960s. As drift term = — 0.44Ax and as diffusion 
term = 0.003exp(-r/2) + 0.019(Ax + 0.04) 2 were taken. The units of Ax 
are multiples of the standard deviation a of Ax with At = 40960s. The pdjs are 
shifted in vertical directions for convenience of presentation. 

Figure 2: Joint pdf p(Ax 2) Ai 2 ; Axi, At\) for the simultaneous occurrence of 
price differences A#i(A£i) and Ax 2 (At 2 ) for US dollar-Deutsch mark exchange 
rates. The contour plot is shown for At\ = 6168s and At 2 = 5120s. a = 
0.0635L>M denotes the standard deviation of x(t). The contour lines correspond 
to \ogp = — 1, . . . , —4. If the two price changes were statistically independent 
the joint pdf would factorize into a product of two pdfs: p(Ax 2) Ai 2 ; Axi, At\) = 
p(Ax 2) At 2 )p{Axi ) Aii). The tilted form of the joint pdf provides evidence that 
such a factorization does not appear for small values of | log(Ati/Ai 2 )|- 

Figure 3 : Contour plot of the conditional pdfp(Ax 2) At 2 |Axi, At\) =p(Ax 2 , At 2 ; Axi, Ati)/ 
p(Axi, Ati) for At 2 = Ati/(1.2) 2 and t x = 5120, in the range (-3a < Ax { < 3a, 
i = 1,2). a denotes the standard deviation of x(i), see fig. 2. The contour lines 
correspond to logp = — 1, . . . , —4. For statistically unrelated quantities the con- 
ditional pdf would reduce to the probability density p(Ax 2) At 2 ). 
In order to verify the Chapman-Kolmogorov equation, the directly evaluated pdf 
(black lines) is compared with the pdf calculated by p ca i{Ax 2) At 2 \Axi, Ati) = 
) dxip{Ax 2) At 2 \Axi, Ati)p(Axi, At i \Ax 1 , Ai x ) for At 2 = Aii/(1.2) 2 ; At* = Ati/1.2 
(red lines). Assuming an statistical error of the square root of the number of 
events of each bin we find that both pdfs are statistically identical. 
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